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The modelling of fluid particle accelerations in homogeneous, isotropic turbulence in 
terms of second-order stochastic models for the Lagrangian velocity is considered. The 
basis for the Reynolds model (A. M. Reynolds, Phys. Rev. Lett. 91(8), 084503 (2003)) 
is reviewed and examined by reference to DNS data. In particular, we show DNS data 
that support stochastic modelling of the logarithm of pseudo-dissipation as an Ornstein- 
Uhlcnbeck process (Pope and Chen 1990) and reveal non-Gaussianity of the conditional 
acceleration PDF. The DNS data are used to construct a simple stochastic model that is 
exactly consistent with Gaussian velocity and conditionally cubic-Gaussian acceleration 
statistics. This model captures the effects of intermittency of dissipation on accelera- 
tion and the conditional dependence of acceleration on pseudo-dissipation (which differs 
from that predicted by the refined Kolmogorov (1962) hypotheses). Non-Gaussianity of 
the conditional acceleration PDF is accounted for in terms of model nonlinearity. The 
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diffusion coefficient for the new model is chosen based on DNS data for conditional two- 
time velocity statistics. The resulting model predictions for conditional and unconditional 
velocity statistics and timescales are shown to be in good agreement with DNS data. 



1. Introduction 



Of late, the statistics of f luid particle accel e ration in turbulence have been the subjec t 



of many experimental (e 



Christensen k, Adrian 



g..lVoth et al 



1998) 



200 21): iGvlfason et a 



La Porta 



2004); 



et al. 



2001) 



Voth et al 



Mordant et al 



(2002); 



(2004 )) and nu- 



merica l (e.g., 



2005) 



2001 



Yejmgl 11997) 



Vedula fc Yeung 



(1999) 



Biferale et al. 



(2005) 



effor t s. These investig ations have spurred a renewed interest l|Pone 



2002); 



Reynold; 



Ygung. 



2002) 



sLuL 



Beck 



12003)) in the modelling of conditional and unconditional statis- 
tics of velocity and acceleration in terms of second-order Lagrangian stochastic models. 
Recent work in this area has focused on the construction of stochastic models that are ca- 
pable of reproducing intermittency and Reynolds number effects in Lagrangian statistics 
as observed in experiments and DNS. In other words, second-order stochastic models 
can be formulated in such a way as to incorporate accurate one-time statistics (and 
their Reynolds-number dependence) from experiments or DNS and be able to reproduce 
intermittent two-time statistics in good agreement with experiments or DNS. 



Reynolds- number effects in Lagrangian stochastic models were first addressed by 



Sawforc 



11991). The Sawford (1991) model is exactly consistent with a joint-normal stationary 
one-time distribution for Z = [U, A] T , where U(t) and A(t) denote (modelled) stochastic 
processes for one component of the Lagrangian velocity and acceleration. As a result, the 
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SDEs for the Sawford (1991) model are linear: 



dZ 






1 




b 2 




2^1 








Zdt + 
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(1.1) 



where ov and cr^ denote standard deviations for velocity and acceleration, b is a dif- 
fusion coefficient, and W is a standard Brownian motion (or Wiener process). Sawford 
(1991) showed that matching of the second-order Lagrangian velocity structure function 
Du(s) — ({U(t + s) — U(t)) 2 ) with the Kolmogorov (1941) hypotheses for the universal 
equilibrium range uniquely identifies the diffusion coefficient as 



2a^(T-- 1 +^ 1 )T--V, 



(1.2) 



where T£° 



Co (e) anC1 tr > 



— '— Here, Co is the Kolmogorov constant for the 



2a y <e> 

second-order Lagrangian velocity structure function, ao is the acceleration variance nor- 
malized by the Kolmogorov scales, (e) is the mean dissipation and v is the kinematic vis- 
cosity. Sawford (1991) also showed that model predictions are very close to DNS data for 
unconditional velocity and acceleration autocorrelations at low Reynolds number. How- 
ever, the Sawford model ignores intermittency of Lagrangian statistics and incorporates 
a Gaussian Lagrangian acceleration P DF, at variance with th e observed n on-Gaussianity 



of acc e leration found in ex periments (|La Porta et 



(1989) 



Ycuneei 



|200lJ)) and DNS JYeung fc Pope 



(2005)). 

Reyno lds ( 2003j) addressed the problem of incorporating a strongly non-Gaussian PDF 
of acceleration into a Lagrangian stochastic model. Reynolds showed that an improved 
representation for the Lagrangian acceleration PDF in a second-order stochastic model 
can be obtained by explicitly accounting for intermittency of dissipation. Specifically, he 
assumed a log-normal distribution for the dissipation rate, e, together with a Gaussian 
assumption for the conditional PDF of A\e. The latter assumption may be restated in 
terms of the conditionally standardized acceleration defined by A = (which has 



4 A. G. Lamorgese 1 , S. B. Pope 1 , P. K. Yeung 2 and B. L. Sawford 3 

zero and unit values for its conditional mean and variance). In the Reynolds model, 

the conditional distribution A\e is assumed to be universal and, in particular, standard 

normal. This may be interpreted to imply that intermittency of dissipation is solely 

responsible for intermittency in acceleration. 

Reynolds also assumed Gaussian velocity statistics and independence of velocity from 
dissipation and acceleration. In other words, the Reynolds model is (by construction) 
exactly consistent with a joint-normal stationary one-time distribution of (U, A, hie). 

To completely specify his model, Reynolds assumed the Kolmogorov (1962) prediction 
for the conditional acceleration variance, 



<J 2 A Jal = al(e/{e)ff\ 



(1.3) 



where clq is a Ko lmogorov cons t ant an d a n = ((e) 3 /^) 1 / 4 is the Kolmogorov acceleration 



scale. Following iPooe fc Chenl l|199(tl . Reynolds also assumed an Ornstein-Uhlenbeck 
(OU) process for \ = me - The resulting model can be written as an SDE for Z = 
[[/,1,1ns- (lne)] T : 



b/<y A \e 



dZ 






°~A\e 
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dW 
dW 



(1.4) 



In these equations, o~ x and T x denote the standard deviation and the integral scale for 
X, whereas W and W are independent Wiener processes. The dissipation equation is 
effectively decoupled from the rest of the system and therefore the Reynolds model is 
linear in U and A. Additional assumptions made by Reynolds are: (i) a choice of diffusion 
coefficient made by analogy with the Sawford (1991) model, i.e., 



b=J^ 2 U (T- 1 £ + t^)T-^ s , 



(1.5) 
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where Tl i6 = ~§^^f~ an d tr/,e — (Co being a model constant), and (ii) T x (e)/afj = 

const. 

In this paper, we first review the basis for the Reynolds model against DNS. Then, a 
novel stochastic model is constructed that incorporates one-time information from DNS 
and yields model predictions for two-time velocity statistics in good agreeement with 
DNS. 

The plan for this paper is as follows. In Section[21we review DNS data (first presented in 



Yeung et al 



(2005)) for intermittency of dissipation, the PDF of conditionally standard- 
ized acceleration and the variance of acceleration conditioned on the pseudo-dissipation. 
In Section|21 a novel stochastic model that incorporates non-Gaussian one-time statistics 
from DNS is formulated. Non-Gaussianity of the conditionally standardized acceleration 
PDF is accounted for in terms of nonlinearity in the model. In Section^ we show a choice 
of diffusion coefficient based on DNS data for conditional velocity autocorrelations that 
yields model predictions for conditional and unconditional velocity autocorrelations and 
timescales in good agreement with DNS. Conclusions for this work are summarized in 
Section |5] 



2. DNS Data for Stochastic Modelling 

2.1. Intermittency of Dissipation 

Figure n shows (one-time) PDFs for me, ln£ and ln</? for R\ « 680 (R\ — y 
being the Taylor-scale Reynolds number), where e = 2vsijSij is the dissipation rate, 
£ = IvTijTij is the "enstrophy" , and ip — vuijmj is the pseudo-dissipation (sij and 
Tij being the strain-rate and rotation-rate tensors, i.e., u^.j = Sjj + r»j). This figure 
suggests that pseudo-dissipation (as opposed to the dissipation rate, or the enstrophy) 
is closest to log-normal for R\ « 680. In fact, DNS data support this conclusion for 




Reynolds numbers in the range R\ sw 140 — 680. In this paper we do not purport to 
discuss the validity of the log-normal model as opposed to more accurate intermittency 
models for dissipation. We limit ourselves to the observation that pseudo-dissipation may 
be approximately described as log-normal for R\ w 140 — 680. 
Let us now assume that x = ^ ni f/( £ ) evolves by an OU process, 

/ a 2 \ dt /24 
dX = -{ X+ f)T x + ^ dW ' ™ 

It follows that (X(t)\X(0) = 4>) = ^e"*/^, where X = (Imp - (ln(p))/a lnlp . Figure^ 

shows (Lagrangian) two-time conditional means of X from DNS (Fig. [21 shows two-time 

conditional means from DNS with X re-defined in terms of e). These figures suggest 

that two-time conditional means of the (standardized) logarithm of pseudo-dissipation 

are closest to simple exponentials (at least for R\ « 140 — 680, based on DNS at different 
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r/T x 

Figure 2. DNS data for conditional expectations (solid) with X = (Ins — (lne))/o"i n 
vs. (/>e~ T / Tx (dashed, Tx being the Lagrangian integral timescale for X). 



Rx ~ 400 




r/Tx 

Figure 3. DNS data for conditional expectations (solid) with X — (In tp — (In ip)) / a\ n 
vs. (f>e~ T ^ Tx (dashed, Tx being the Lagrangian integral timescale for X). 
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Rx 38 140 240 400 680 



T x /T 0.3317 0.1719 0.1311 0.1056 0.0906 



Table 1. Values of T x /T (with T = 1.5<^/<e)) from DNS of homo geneous turbulence. 



Reynolds numbers). Therefore, on the basis of DNS, we argue that, for the purposes of 
stochastic modelling, pseudo-dissipation is well-approximated by an OU process. 

To completely specify the dissipation model 12.1]) (wher e \ = m( ^/( e ))j °~x anc ^ -^x 



have to be prescribed. As discussed bv lYeung et all (|2005f> . the DNS data support the 
Kolmogorov (1962) prediction for 



a 2 x = A + ^-hxRx, 



(2.2) 



with /i « 0. 25, in good agreeme nt with 



reported by 



Yeung et al 



Sreenivasan fc Kailasn ath (1993) (A ps —0.863 is 



(2005J)). The integral timescale T x is chosen to match the DNS 



data (Table P). 

2.2. PDF of Conditionally Standardized Acceleration 
We now investigate the conditionally Gaussian assumption in the Reynolds model. Figure 
4 shows PDFs f^ e (a\e) for different values of i for Rx ~ 680 (a standard normal PDF 
and the unconditional acceleration PDF for R\ « 680 are also shown). Figure 5 shows 
analogous information when cp is used in place of e. Both sets of conditional PDFs are 
much less intermittent (with weaker tails at large fluctuations) than the unconditional 
acceleration PDF. This is particularly true of the PDFs of A\ip (with A = A/aA\ v ) which, 
however, still show significant non-Gaussian behaviour. A remarkable degree of collapse 
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-60 -45 -30 -15 15 30 45 60 

a/<?A\e, a/ a a 

Figure 4. Standardized onditional PDFs of A given e = i from 2048 3 DNS at 
Rx ~ 680. Lines A-E are for e/(e) = {0.0359, 0.136, 0.469, 1.62, 6.05}, corresponding to 
Z = (lne- (e))/crine = {-2.054, -0.994, 0, 0.994, 2.054}. Also shown are the unconditonal 
PDF of acceleration (solid unmarked line) and a standard Gaussian PDF (dashed). 

of these PDFs for different values of the conditioning variable is notable (except for very 
small and very large conditional fluctuations). Therefore, to a first approximation, the 
PDFs /^(alv?) may be described as approximately independent of (p. In fact, based on 
simulations at different Reynolds numbers (not shown), the conditional PDF of A\ip may 
also be (approximately) described as independent of the Reynolds number. 
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Q 




15 15 



45 60 



Figure 5. Conditional PDFs of A given ip = <p from 2048 3 DNS at R x w 680. 
Lines A-E are for ip/((p) = {0.0362, 0.134, 0.458, 1.56, 5.79}, corresponding to 
Z = (\nip - (ip})/cr lnv = {-2.054, -0.994, 0, 0.994, 2.054}. Also shown are the uncondi- 
tonal PDF of acceleration (solid unmarked line) and a standard Gaussian PDF (dashed) . 



Yeung et all 1)20051) suggest that the PDF of A\ip can be described (to a very good 
approximation) as cubic-Gaussian. By definition, a random variable Z is cubic- Gaussian 
with parameter p (also denoted as Z ~ G 3 (p)) if 



Z = C[{l-p)X+pX\ 



(2.3) 



where X is a standardized Gaussian random variable and C is determined by the stan- 
dardization condition (Z 2 ) — 1 as C(p) = (1 + 4p + lOp 2 )^ 1 / 2 . Figure HJ1 shows that the 
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Figure 6. Standardized conditional PDFs of A given 9? = (p. The symbols are 
from the 2048 3 DNS; the lines are the cubic-Gaussian PDF with the same kur- 
tosis as the DNS data. The values of the conditioning variable are such that 
Z = (\n<p - {(p))/a ltllp = {-2.054, -0.994, 0, 0.994, 2.054}. In each plot, the lower curve and 
the j/-axis correspond to the lowest conditioning value. The curves for the other conditioning 
values are successively shifted upwards, by an amount 0.2 in the linear plot (left), and by a 
factor of 100 in the logarithmic plot (right). 



cubic-Gaussian PDF provides a remarkably accurate description of the conditional PDFs 
f^tpi^W)- Comparable accuracy to that in Fig. is achieved when the cubic-Gaussian 
PDF is use d to fit DNS data for fx^i&W) a t lower Reynolds numbers (not shown, but see 



Sect. IV of 



Yeung et al. 



l2005j)). A value of p ~ 0.1 results from the observation (based 



on DNS) th a t U4(A \(p) ss 8 (approximately independent of ip and R\; more details in 



Yeung et al. 



(2005)). 
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Figure 7. Variance of acceleration conditioned on the pseudo-dissipation for different values of 
R\. The symbols are the DNS data; the lines are the empirical fit Eq. 12.511 . The lowest curve 
and the y-axis correspond to R\ m 38. The other four curves are for R\ m 140, 240, 400, 680, 
successively shifted upwards by a factor of a/10- 



2.3. Conditional Acceleration Variance 



With a view to the joint (stochastic) modelling of acceleration and pseudo-dissipation, 
we now investigate the validity of the Kolmogorov (1962) prediction for the conditional 
acceleration variance, 

a 2 Alv /a^a*^/(e)f/ 2 . (2.4) 

Figure shows ^-dependences from DNS for the conditional acceleration variance at 
differe nt Reyno l ds nu mbers. In the same figure, the following expression (first presented 



Yeung et al 



(2005)) 



20 



(2.5) 
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is shown to be an accurate representation (except at the smallest R\) of the DNS data. 

As may be seen, the low- tp behaviour for the conditional acceleration variance deviates 

strongly from that predicted by Eq. I|2.4() . Also, careful measurement of the slope for the 

large-(p portion of the curves in Fig.[7|yields values that are systematically less than 1.5, 

again at variance with the Kolmogorov (1962) prediction. 

Equation 1)2. 5|1 is most useful for stochastic modelling purposes because it accurately 

parameterizes the conditional acceleration variance (given ip) in terms of both the value 

being conditioned upon and the Reynolds number. 

3. Conditionally Cubic-Gaussian (CCG) Stochastic Lagrangian 
Models 

Lagrangian statistics for e, £ and tp from DNS show that stochastic modelling is easiest 
when pseudo-dissipation is used in place of the dissipation rate or the enstrophy. This 
is because (i) tp is closest to log-normal, and (ii) two-time conditional means of In tp are 
closest to exponential, and (iii) the conditional PDFs of acceleration given tp = tp collapse 
best, and with the narrowest tails. Thus, we base the model on pseudo-dissipation tp and 
take the OU process Eq. 1)2. lfl as its stochastic model. 

Conditioning on pseudo-dissipation is most useful when considering the joint-statistics 
of acceleration and pseudo-dissipation because the PDF of A\tp may be described (to 
a first approximation) as universal and, in particular, cubic-Gaussian. In other words, 
given ip = tp and a standardized Gaussian random variable A, the acceleration A can be 
modelled as 

A = a Al0 C[(l~p)A+pA 3 }. (3.1) 

For given tp, the relation between A and A is invertible (and one-to-one). 

The stochastic model is most conveniently expressed in terms of the velocity U(t), 
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the "Gaussian" acceleration A(i) (related to the acceleration by Eq. (f 3 . 1 1> ^ - and the log- 
pseudo-dissipation variable x* = X ~ (x) = ( e )) — {^ n ( l fi/ ( e )))- 
The model is 

dU = Adt = a Alv C[(l - p)A + pA 3 ]dt, (3.2) 
dA = M + bdW, (3.3) 

where 9 and b are drift and diffusion coefficients specified below. 

The stationary one-time joint PDF of U, A and x* IS denoted by f(v, a,x*), where v, a 
and x* are sample-space variables corresponding to U, A and x* ■ We now assume that 
this PDF is joint-normal with the variables being uncorrelated with each other (at the 
same time). Thus, with the assumptions made the joint PDF is 

f = ^ exp vfc exp ^ exp irh) ■ (3 - 5) 

The imposition of this PDF leads to a constraint for the drift coefficient 9 in Eq. (|3.3(l . 
namely, 

9(v,d, V ) = -^Cv(l+p + pd 2 ) + -—lTib 2 f + 9*, (3.6) 
afj 2 aa 

where 9* is any function such that -§g{9* f) = 0, which for simplicity we take to be zero. 
We also introduce the assumption db/dd — (which can be supported using an adiabatic 
elimination argument in the limit R\ — > oo). Then, Eq. 1)3.3(1 can be rewritten as 

dA = —Adt- '^UC(l + p + P A 2 )dt + bdW. (3.7) 
2 a v 

This equation (together with Eqs. ((3.2(1 and ((3.4(1 ^ defines a class of CCG models, 
i.e., different models with the same stationary distribution 1(3. 5(1 correspond to differ- 
ent choices of b. Each model captures the conditional dependence of acceleration on 
pseudo-dissipation based on DNS (Eq. 1(2.5(1 ) that accounts for deviations from the Kol- 



Conditionally Cubic- Gaussian Stochastic Lagrangian Model 15 
Rx w 680 

1 

0.9 
0.8 
0.7 

8 0-5 

P-H 

o 

8 0.4 

-(^ 

t 0.3 

I 0.2 

> 

0.1 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

t/T 

Figure 8. Velocity autocorrelation for CCG model (solid) vs. DNS data (dashed). The 
thicker, dashed line is for the Reynolds (2003) model. 

mogorov (1962) hypotheses. Also, each model is non-linear because it accounts for the 
non-Gaussianity of the conditionally standardized acceleration PDF. 

4. Specification of Diffusion Coefficient 

The diffusion coefficient for the Reynolds (2003) model is specified by analogy with 
the Sawford (1991) model (Eq. p. 5(1 ). This choice is arbitrary and, furthermore, it leads 
to poor agreement with DNS for unconditional autocorrelations (Fig. [BJ. 

We now investigate the question of how to select b(x) by reference to two-time condi- 
tional velocity statistics from DNS. To this end, the Reynolds model with a frozen dissi- 
pation x = X is considered. Then, the first two equations in (|1.4f> revert to a Sawford '91 
type of model. With the frozen model, we investigate x-dependences for b such that the 
frozen model matches the DNS values of o~a\ x and t\j\ x (ju\ x being halving times for con- 
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R. = 240 

A 




Figure 9. Phi-dependences for diffusion coefficient (solid) for the Reynolds (2003) model with a 
frozen dissipation that matches the DNS values of <?a\x and Txj\ x . The dashed line is a tentative 

/ \ 0.3 

fit to the solid curves G(ip) = ( ) (Eq. Kll ). 



ditional velocity autocorrelations pu\ x (t) = (U (t)U (0)\x{0) — x): i- e -j Pu\x( T u\x) = V^)- 
Results from such calculations for different values of R\ are shown in Fig. This figure 
would suggest that (^-dependences for (r^b 2 / a 2 ) / {o~\\ x / a 3 ) at different Reynolds num- 
bers may be approximately described in terms of a single function G(ip) = I 7jy J • The 
diffusion coefficient is then given by 




where u v — (v(e)) 1 ^ 4 is the Kolmogorov velocity scale. We then use this determination 
to specify b in (|3.3|) . i.e. we let 

b 2 = c(R x )G^)^, (4.2) 

Ur, 

where c(R\) is a correction factor (given in Table[2l so selected as to ensure good agree- 
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Figure 10. Conditional velocity autocorrelations for CCG model (solid) vs. DNS (dashed). 



ment between model predictions and DNS data for unconditional velocity autocorrela- 
tions (Fig. [HJ. The resulting model predictions for conditional velocity autocorrelations 
and timescales are shown in Figs. ^| and ^] Model predictions for unconditional and 
conditional acceleration autocorrelations are shown in Figs. ^| and These plots sug- 
gest that the specification Ij4.2|l is capable of reproducing reasonable agreement with 
DNS data for conditional and unconditional velocity autocorrelations and timescales. 
Also, Eq. 14.2fl yields PDFs of Lagrangian velocity increments that are approximately 
Gaussian for large time-lags (Fig. ll4|) . These PDFs develop stretched tails as the time-lag 
decreases and ultimately approach the Lagrangian acceleration PDF for very small time- 
lags. This behaviour is consiste nt with recent observa tions of Lagrangian intermittency 



in experiments and simulations (Mordant et al. 



(2004)). 
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Figure 11. Model predictions for halving times for conditional velocity autocorrelations 

(solid) vs. DNS (dashed). 



Rx 38 140 240 400 680 



c(R x ) 1.4 1.0 1.0 1.0 0.7 



Table 2. Values for c(R x ) in Eq. tOl . 



5. Conclusions 

After a brief review of the basis for the Reynolds (2003) model against DNS, we 
have shown the formulation of a novel stochastic based on simple data assimilation from 
DNS. This model is exactly consistent with Gaussian velocity and conditionally cubic- 
Gaussian acceleration statistics and incorporates a representation for the logarithm of 
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Figure 12. Acceleration autocorrelation for CCG model (solid) vs. DNS (dashed). The 
thicker, dashed line is for the Reynolds (2003) model. 
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Figure 13. Conditional acceleration autocorrelations for CCG model (solid) vs. DNS (dashed). 
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( o i R x g 680 




v/au, cl/oa, etc. 



Figure 14. Model predictions for standardized PDFs of Lagrangian velocity increments. 



pseudo-dissipation as an Ornstein-Uhlenbeck process. The new model captures the effects 
of intermittency of dissipation on acceleration and the deviations from the Kolmogorov 
(1962) hypotheses (based on DNS) in the conditional dependence of acceleration on 
pseudo-dissipation. Further, non-Gaussianity of the conditionally standardized accelera- 
tion PDF (as observed in DNS) is captured in terms of model nonlincarity. An empirical 
specification of diffusion coefficient based on DNS data for conditional velocity timescales 
yields reasonable agreement with DNS data for conditional and unconditional velocity 
autocorrelations and timescales. 
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